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A bird distribution model for ring recovery data: where do 
the European robins go? 
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For the study of migratory connectivity, birds have been individually marked 
by metal rings for more than 100 years. The resulting ring recovery data have 
been compiled in numerous bird migration atlases. However, estimation of 
what proportion of a particular population is migrating to which region is con- 
founded by spatial heterogeneity in ring recovery probability. We present a 
product multinomial model that enables quantifying the continent-wide distri- 
bution of different bird populations during different seasons based on ring 
recovery data while accounting for spatial heterogeneity of ring recovery proba- 
bility. We applied the model to an example data set of the European robin 
Erithacus rubecula. We assumed that ring recovery probability was equal 
between different groups of birds and that survival probability was constant. 
Simulated data indicate that violation of the assumption of constant survival 
did not affect our estimated bird distribution parameters but biased the esti- 
mates for recovery probability. Posterior predictive model checking indicated a 
good general model fit but also revealed lack of fit for a few groups of birds. 
This lack of fit may be due to between-group differences in the spatial distribu- 
tion on smaller scales within regions. We found that 48% of the Scandinavian 
robins, but only 31% of the central European robins, wintered in northern 
Africa. The remaining parts of both populations wintered in southern and cen- 
tral Europe. Therefore, a substantial part of the Scandinavian population 
appears to leap over individuals from the central European population during 
migration. The model is applied to summary tables of numbers of ringed and 
recovered birds. This allows us to handle very large data sets as, for example, 
those presented in bird migration atlases. 
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Introduction 

The study of migratory connectivity is important to 
understand the ecology and population dynamics of spe- 
cies (Webster et al. 2002). Ecological conditions encoun- 
tered during the nonbreeding period can affect breeding 
performance due to carry-over effects (Norris and Marra 
2007; Schaub et al. 2011; Harrison et al. 2013). For con- 
servation of bird species, it is therefore valuable to know 
where individual birds are found throughout the course 
of a year. In addition, knowledge of migratory connectiv- 
ity is important to assess the risk of avian transmitted 
diseases such as avian flu (Liu et al. 2005). 

There are many challenges to studying the migratory 
pathways of birds, particularly smaller species. An ideal 



technique to study bird migration would be a device that 
is so small that it could be carried by small migratory 
species without affecting their behavior, while allowing 
long-distance tracking with high precision (Robinson 
et al. 2010). However, such a device does not yet exist 
(Robinson et al. 2010). The existing tracking techniques 
are subject to a trade-off between weight and temporal 
and spatial precision of the measurements. 

GPS devices measure location with high precision (a 
few meters, e.g., Hiinerbein et al. 2000). A GPS can be 
combined with a satellite transmitter to follow the track 
of a bird in real time (e.g., Kjellen et al. 1997). However, 
the high weight of GPS and satellite telemetry devices 
(normally >10 g, Robinson et al. 2010) prohibits the use 
of these techniques for birds smaller than around 200 g 
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body mass (usually an additional weight of <5% of the 
body mass is accepted, e.g., Naef-Daenzer et al. 2001). 
Even for larger species, researchers often only equip those 
individuals that are heavier than a threshold weight (e.g., 
Sauter et al. 2012). This produces a nonrandom sample 
of the population being studied. They are also very 
expensive, limiting the number that can be deployed. 

Lightweight-tracking devices are data loggers such as 
light-based geolocators (e.g., Stutchbury et al. 2009; 
Bachler et al. 2010). The smallest geolocator now weighs 
around 0.6 g (Bridge et al. 2013) and can be applied to 
birds with a body weight of at least 10 g. However, the 
derivation of the location from light measurements is 
associated with large errors (up to several 100 km) 
depending on the season, the behavior of the bird and 
the shading the bird has experienced during its journey 
(Lisovski et al. 2012). Furthermore, as the birds have to 
be recaptured in order to retrieve the data, only those 
individuals that survive and return to the place of capture 
(normally the breeding site) can be tracked. 

As an alternative to mounting a tracking device on the 
bird's body, analyzing the chemical composition of 
the feather can give information on the whereabouts of 
the bird during the time when the feathers grew. Com- 
monly, isotopes (Hobson and Wassenaar 2008) or trace 
elements (Szep et al. 2003) are used for this purpose. 
However, the spatial resolution is normally very low and 
often it is only possible to determine whether or not two 
individuals used similar wintering areas. 

In contrast, the most widely used method of individu- 
ally marking birds, using metal rings, can be applied to 
nearly all bird species, and it provides precise location 
information if the bird is found and its ring is reported 
to a ringing scheme. Marking birds with metal rings has 
been performed by many professionals and amateur orni- 
thologists for more than a century (Mortensen 1901). 
Since that time, more than 100 million birds have been 
ringed, and several million recovery records have been 
collected in international databases such as the Euring 
database (www.euring.org). Given the large sample sizes 
and wide geographic area of coverage, these data would 
seem to be ideal for studying migratory connectivity, 
defined as the proportions of birds from different breed- 
ing populations migrating to different wintering areas 
(Webster et al. 2002). 

However, ringing data present challenges for analyzing 
migratory connectivity, because the probability of finding 
and reporting a ring varies geographically and over time, 
resulting in a nonrandom sample of the population being 
studied (Perdeck 1977). 

Numerous bird migration atlases (Zink and Baiiiein 
1995; Brewer et al. 2000; Fransson and Pettersson 2001; 
Wernham et al. 2002; Bakken et al. 2003; Bonlokke et al. 



2006; Spina and Volponi 2008; Csorgo et al. 2009; 
Hiippop and Hiippop 2009; Saurola et al. 2013) have 
mapped data that link ringing with recovery locations. 
These bird migration atlases show where each bird species 
can go, but they do not allow us to quantify migratory 
connectivity, because the probability of finding a ringed 
bird and reporting its ring to a ringing scheme varies tre- 
mendously among different regions (Perdeck 1977; Fiedler 
et al. 2004; Clark et al. 2009; Korner-Nievergelt et al. 
2010a). None of the migration atlases we know of have 
made an attempt to formally take into account recovery 
probability in their recovery maps. 

A few approaches have been proposed for estimating 
geographic variation in recovery probabilities, using simu- 
lation (Lokki and Saurola 2004), using covariates for 
recovery probability (Cowen this volume) or by compar- 
ing groups of birds for which equal recovery probability 
is assumed (Busse and Kania 1977; Jenni 1987; Kania and 
Busse 1987). The latter method has been used in recovery 
models that aimed to describe migratory connectivity 
(Bauthian et al. 2007; Thorup and Conn 2009; Korner- 
Nievergelt et al. 2012). 

Here, we present a further alternative of a large-scale 
spatial ring recovery model. The aim is to describe a 
model formulation that can potentially be used to 
account for spatially and temporally heterogeneous ring 
recovery probability in ring recovery maps for future 
migration atlases. Therefore, the model is simple so that 
it can potentially be applied to a large number of data 
sets from various species. 

We apply the model to example data from the Euro- 
pean robin Erithacus rubecula. The European robin is a 
partial migrant. It breeds all over Europe. The Scandina- 
vian population migrates southwards. In southern and 
central Europe, the wintering area overlaps with the 
breeding area. The southernmost wintering sites are in 
northern Africa. However, little is known about what pro- 
portion of each population migrates and how far the 
migrants go. 

Methods 

Data 

We used ringing and recovery data for the European 
robin, which were kindly provided by the ringing stations 
Ottenby (www.sofnet.org/ottenbyfagelstation/start), Falst- 
erbo (www.falsterbofagelstation.se), and Christianso 
(Denmark) and the ringing schemes Hiddensee (Germany) 
and Switzerland (Fig. 1). Only data from birds ringed 
when fully grown were used. We used only reencounters 
of dead birds (i.e., recoveries) found more than 5 km 
from the place of ringing. This selection criterion ensures 
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Figure 1. The four regions defined in this study. The circles identify 
the ringing locations. 

that the proportion of sedentary birds is not overesti- 
mated because ring recovery probability can be enhanced 
close to the ringing places (due to the activity of the ring- 
ers or the awareness of the people). Furthermore, we 
selected those recoveries for which the finding date was 
known with a precision of a least ±2 weeks. Recoveries 
from the month of ringing were discarded (i.e., only the 
first encounter per month was used). We only used data 
for birds ringed in 1964 or later, to ensure that all data 
sets were from a similar time period (Table 1). 



We divided Europe into four regions: (A) Fennoscandia 
including Chistianso, which is an island in the Baltic Sea 
south of Sweden; in this region, the European robin is a 
breeding or migratory bird, but it does not overwinter; 
(B) central Europe, which is both a breeding and a win- 
tering area; (C) southern Europe; and (D) northern Africa 
(Fig. 1). Southern Europe (C) and northern Arica (D) 
can be considered as purely wintering areas for the Euro- 
pean robins ringed in northern and central Europe. The 
robin also breeds in southern Europe and along the 
northern edge of northern Africa (Bauer et al. 2005). 
However, we assumed that we have no such birds in our 
data as all birds were ringed north of these regions. 

We defined 12 release occasions according to the 
month of release. For each month, we summed data 
across years. Of the birds released in Fennoscandia, we 
can assume that they bred in Fennoscandia or in the 
northwestern part of Siberia. Among birds released in 
central Europe, only those released between May and July 
can be considered as breeding in central Europe. Birds 
ringed in central Europe during autumn, winter, and 
spring (August to April) have to be considered as a mix- 
ture of central European and Scandinavian birds as the 
Scandinavian birds migrate through or winter in central 
Europe during these months (Bauer et al. 2005). 

Model 

We used a multistate model as described by Arnason (1972) 
and Schwarz (1993). We had four states representing the 
regions (k = A, B, C, D). The model described by Arnason 
(1972) and Schwarz (1993) estimates the transition 



Table 1. Description of the five data sets analyzed in this study with the total numbers of ringed European robins, the numbers of ringed birds in 
each month summed over the years, and the total number of ring recoveries. 



Place of ringing 


Ottenby (S) 


Falsterbo (S) 


Christianso (DK) 


Hiddensee (D) 


Switzerland (CH) 


Region 


Fennoscandia 


Fennoscandia 


Fennoscandia 


Central Europe 


Central Europe 


Years 


1964-2011 


1980-2011 


1 994-2002 


1964-2011 


2008-2011 


Total numbers ringed 


200,230 


111,421 


28,006 


187,693 


17,596 


Numbers ringed in: 












January 


0 


0 


0 


222 


8 


February 


0 


0 


0 


173 


2 


March 


2100 


1650 


135 


7704 


590 


April 


63,206 


26,890 


8961 


40,328 


2099 


Mai 


20,163 


5227 


1441 


5672 


6 


June 


142 


52 


37 


2888 


0 


July 


42 


14 


0 


5690 


28 


August 


B20 


2410 


149 


8258 


436 


September 


54,348 


46,072 


10,255 


65,214 


5466 


October 


58,069 


28,441 


6789 


48,729 


8671 


November 


1633 


665 


239 


2418 


270 


December 


7 


0 


0 


397 


20 


Total number of recoveries 


366 


185 


45 


540 


16 
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probability that an individual moves to any possible state k 
at occasion q given it was in a specific state k at occasion 
q — \{m kqlkq ) . However, we formulated the model so that 
these transition parameters were independent of the state of 
an individual at occasion q-l but instead depended on the 
group (or set) to which the bird belonged. Thus, in our 
model, these parameters represented the proportion of a set 
of birds (released in region i during month ;') that were in 
region k during season q (tnjj kq , see Table 2 for notation). 
We assigned the birds to 24 sets based on the ringing loca- 
tion (region A or B) and months (see Table 1 for number 
of ringed birds). The distributions of these sets of birds 
among the four regions were estimated for 8 seasons: winter 
(December to February), March, April, May, summer (June 
to August), September, October, and November. We 
assumed that the robins do not migrate during the three 
winter months and the three summer months, respectively. 
Note that for grouping the birds released, we use months 
here (indexed by;') but the distribution parameters are esti- 
mated per season (recovery occasion, which is sometimes a 
month and sometimes a group of months, indexed by q). 

We formulated one multinomial model for each set of 
birds ringed in the same region and during the same 
month ;'. 



Table 2. Indices, parameters, and notation for data. 
Data 

R,j Vector of length K*Q + 1 containing the number of 

recoveries in each region (k) and season (q) of the set of 
birds released in region / during month /'. The last number 
in the vector is the number of birds never recovered. 

N,j Number of ringed and released birds at region / during 
month j. 

I Total number of release regions (2) 
J Total number or release occasions, that is, months (12) 
K Total number of recovery regions (4) 
0 Total number of recovery occasions, that is, seasons (8) 
Indices 

/ Release region, Fennoscandia and central Europe (A, B) 
j Release occasion, 12 months (1,..., 12) 
k Recovery region, destination region (A, B, C, D) 
q Recovery seasons: 1....8: winter, March, April, May, summer, 
September, October, November 

f Month (1 12) 

Model parameters 
Pij Probability vector of length K*Q + 1 
Pijk q The first K*Q elements of p,y 

rriijkq The proportion of the set of birds /; (ringed in region / during 
month j) being in region k during season q 

Fj q Probability that a bird alive during month j dies during a 
subsequent season q 

s Monthly survival probability 

r kq Recovery probability; probability that a ringed bird that has 
died in region k during season q is found and reported 



Ry ~ Multinom(p-, Nij) 

This yielded I*J = 2*12 = 24 models (one for each set 
of birds released), which together formed a product 
multinomial model. 

Ry is a vector of length K*Q + 1 containing the num- 
ber of recoveries in each region k and season q for each 
set of birds {if). The last element of the vector contained 
the number of birds that were never recovered. 

The elements of the vector p y and the probabilities pij kq 
are the product of three probabilities: (1) the probability that 
a bird is in region k during season q given it was ringed in 
region i during month ;', m,^, (2) the probability that a bird 
ringed and released during month ;' dies during a subsequent 
season q, F jq , and (3) the probability that a dead ringed bird 
is found and its ring reported to a ringing scheme, r kq . 

Pijkq = tnijkq * Fjq * rkq 

The probability that a ringed bird is never reported; 

K Q 

Pij(K*Q+l) = 1 - ^2^2 Pijkq 
k=l q=l 

is added as the last element to p y . 

The probability that a bird ringed and released during 
month ;' dies during any subsequent season q (of any sub- 
sequent year), Fj q , is obtained from the probabilities that 
a bird ringed and released during month ;' dies during a 
specific month f (in any subsequent year), F*. In our 
model, we ignore the year of recovery and instead only 
consider the month of recovery. If we assume that 
monthly survival s is constant over time and equal 
between the regions, three different cases have to be 
distinguished for the calculation of the probability F? t . 
First, the month of ringing ;' is equal to the month 
of death t. In this case, the bird can survive 11 months 
and die in the 12th or it can survive 23 months and 
die in the 24th, etc. This gives a probability of 
F* = S u (l — s) + S 23 (l — s) + . . .. This is a geometric ser- 
ies that can be formulated as given in the first row of the for- 
mula below. When the bird dies later or earlier in the year 
than it was ringed, this probability has to be adjusted 
according to the difference in the number of months 
between the month of death and the month of ringing. 
These formulas are given in the second and third row below. 



n 



( s"(l-s) 

.G-t-i)(i_ s)+s (/- t )£M 
iS (n-f+;)(l_ s ) +s (i2-m)f^ 



if; = t 
if t > j 
if f <;' 



In order to obtain the probability that a bird ringed during 
month ;' dies during season q (instead of month f), we 
summed the corresponding Fp that is, F } 



tec — ^December 
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"■"-^January ~>~ ■'/February' ana "/summer — -fj Julle + -f/July "+" ■'/August - 

For the recovery seasons of one month duration (March, April, 
May, September, October, November) Fj q = F* . 

The advantage of such a parameterization of survival in 
our model is that we do not need to add year of recovery 
as a fifth dimension (of length 48 in our case) in our 
data. This reduces computing time and enables data to be 
added from additional years without increasing the 
dimension of the data. One drawback is that we lose 
information about temporal variation in survival, but very 
large data sets would be needed to estimate that with 
precision anyway. 

The model was fitted using Bayesian methods. A uni- 
form distribution, Unif(0,l), was used as prior distribu- 
tion for monthly survival probability s. We allowed for 
between-season variance in recovery probability, r kq , with 
partial pooling within each region. To do so, we modeled 
r kq as a realization from a region-specific beta distribu- 
tion: r kq ~ Beta(aj t , fV k ) with Gamma priors for d\ and jV k 
~ Gamma(0.01, 0.01) and P[ ~ Gamma(0.01, 0.01). 

For the bird distribution parameters my kq , we included 
the following constraints. Birds ringed during the three 
winter months (December-February) behaved similarly, 
thus m ilkq = m i2kq = m il2kq , and birds ringed during the 
months June and July were also assumed to behave simi- 
larly, that is, m i6kq = m i7kq . 

We used uniform prior distributions for m^ klj and 
transformed these so that the m ijkq summed to 1 within 
one recovery season for each set of birds, rriijkq = 

+ m hi + m m + w k) ■ For most m U ' we 
used the Unif(0,l) prior distribution. However, because 

most of the European robins have left Sweden by Septem- 
ber and are still in southern Europe in March (Fransson 
and Pettersson 2001), we used Unif(0, 0.01) as prior dis- 
tributions for the proportions of birds in Fennoscandia 
between November and March. In this way, we con- 
strained the proportions of birds wintering in Fennoscan- 
dia to be less than around 1%. Similarly, we constrained 
the proportions of birds spending the summer in north- 
ern Africa and southern Europe during summer to be less 
than around 1% for all sets of birds. 

The model was fitted in JAGS using the package R2jags 
(Su and Yajima 2012). Two Markov chains of length 
20,000 were simulated. Burnin was set to 5000, and the 
chains were thinned by two. Convergence was assessed 
visually and based on the Brooks-Gelman-Rubin statistics 
(Brooks and Gelman 1998). The JAGS code for the model 
is provided in the Data SI. 

In our model, the bird distribution parameters are 
identifiable if there are at least as many sets of released 
birds (if) as there are recovery regions and if the different 
sets of birds differ somewhat in their distribution among 
the regions (Korner-Nievergelt et al. 2010b). Here, we 



have 24 sets of birds and four regions. Thus, all bird dis- 
tribution parameters should theoretically be estimable. 
However, if the data are too sparse or the spatial distribu- 
tion of the different sets of birds is too similar, estimabil- 
ity of some parameters can be poor (Schaub 2009). 
Fitting the model to simulated data sets of different sam- 
ple sizes indicated that for our study, some of the param- 
eters may be little informed by the data (Data S2). 
Therefore, we compared the posterior distribution of each 
parameter with its prior distribution to assess how 
strongly the parameters were informed by the data. For 
this, we calculated the percentage by which the posterior 
distribution overlapped with the prior distribution 
(Garrett and Zeger 2000; Gimenez et al. 2009). Garrett 
and Zeger (2000) suggested that a parameter estimate can 
be considered well informed by the data when the overlap 
between the posterior and a uniform prior is less than 
35%. For nonuniform priors, as used here for the param- 
eters mij kq , the threshold may be higher but no general 
guideline is available (Gimenez et al. 2009). Therefore, we 
used the overlap between the posterior and the prior dis- 
tribution as a relative measure: The smaller the overlap, 
the more data-informed the parameter. 

In order to assess model fit, we used predictive model 
checking (Gelman et al. 1996). To do so, we simulated 
recovery data based on the model while taking into 
account the uncertainty in the parameter estimates. We 
then graphically compared these predictive distributions 
with the corresponding numbers of observed recoveries. 

Results 

Estimated monthly survival probability, s, was 0.89 (95% 
credible interval (CrI): 0.87-0.91, prior overlap: 6%). 
Recovery probabilities, r kq , were lowest in northern 
Africa and highest in central Europe (Table 3). During 
summer, the recovery probabilities were lower than dur- 
ing the migration periods. The mean overlap of the prior 
with the posterior distribution of r kq was 1% (range 
0.1-5%). 

Overlap between prior and posterior distributions for 
the bird distribution parameters m,^ ranged from 27 to 
100% (mean 85%). The estimates were well data- 
informed for Scandinavian birds ringed during April, 
May, September, and October and for the central Euro- 
pean birds ringed between March and October (Fig. 2, 
Data S3). Estimated proportions of birds spending the 
winter in northern Africa were higher for birds ringed in 
September than for birds ringed in October (Fig. 2). The 
average proportion of Scandinavian birds wintering in 
northern Africa is 48% (95% CrI: 36-60%, including only 
estimates that are well informed by the data, i.e., for birds 
ringed during April, May, September, and October 
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Fig. 3). Of the birds ringed in central Europe during the 
breeding season (May to August), the average proportion 
that winters in northern Africa is significantly lower, that 
is, 31% (22-42%). 

Overall, the model seems to fit the data well as assessed 
by posterior predictive model checking (Data S4). Two 
discrepancies between the model and the data deserve 
attention. For Scandinavian birds ringed in October, the 
model predicts a higher than observed number of recov- 
eries in Fennoscandia during summer. We also find more 
recoveries during winter in central Europe of birds ringed 
in central Europe during December than predicted by the 
model. 

Discussion 

The key assumption of the model presented here is that 
the (region (k)- and season (q)-specific) recovery proba- 
bility does not depend on the region (z) and month (f) of 
ringing, that is, the different sets of birds {if) are assumed 
to experience the same recovery probabilities given they 
are in the same region. This assumption could be violated 
if recovery probability is spatially heterogeneous within a 
region and the different sets of birds use different subre- 
gions within a region. For robins ringed in Sweden and at 
the southern Baltic coast, it has been stated that early 
migrants winter more westerly within the Mediterranean 
than late migrants (Hogstedt and Persson 1971; Petters- 
son et al. 1990; Remisiewicz 2002). Furthermore, ring 
recovery probability decreases toward the east in the 
Mediterranean (Korner-Nievergelt et al. 2012). Therefore, 
robins ringed in September (early migrants) will winter 
more to the west and experience a higher recovery proba- 
bility than robins ringed in October (late migrants). 
However, our model assumes the same recovery probabil- 
ity for both sets of birds as we have not divided the Med- 
iterranean into an eastern and western region (Fig. 1). As 
a consequence, our model may underestimate the propor- 
tions of birds ringed in October that winter in southern 
Europe and northern Africa. 

This implies that the difference in the proportion of 
birds wintering in northern Africa between the birds 
ringed in October and those ringed in September may 
appear because the later migrants winter more to the east 
than the earlier migrants (rather than because earlier 
migrants winter farther south). It would be valuable to 
distinguish the eastern and western Mediterranean as sep- 
arate regions in the model. However, when we fitted such 
a model, most of the parameter estimates were not or 
very weakly informed by the data. Therefore, more data 
would be needed to fit such a model. 

Another assumption made by the model is constant 
survival. For two migratory songbirds of similar size to 



) 2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



725 



Bird Distribution Model for Ring Recovery Data 



F. Korner-Nievergelt ef al. 



Released in Fennoscandia 



Released in central Europe 




wmm 



+ 



liiH 




II 1 




23456789 10 12 



Month of release 
Fennoscandia ■ Central Europe 



123456789 10 12 

Month of release 
Southern Europe □ Northern Africa 



Figure 2. Estimated distribution during eight 
seasons (rows) of birds ringed in Fennoscandia 
and central Europe. The month of ringing 
(release) is given on the x-axis. The colors in 
the bars give the distribution of the birds 
among the four regions. Gray bars indicate 
distribution estimates for which the median 
overlap between the prior and the posterior 
distributions was higher than 95%. For these 
groups, the distribution estimates are not 
given. The vertical lines give the 95% credible 
intervals of the summed proportions as given 
in the figure. The 95% credible intervals for 
each single estimated proportion are given in 
the Data S3. For sample sizes, see Table 1 . 



the robin, the black-throated blue warbler Dendroica 
caerulescens and the barn swallow Hirundo rustica, it has 
been shown that survival is reduced in the first few weeks 
after fledging, but later it is similar to that of adults 
(Sillett & Holmes 2002; Griiebler and Naef-Daenzer 2010). 
We excluded birds ringed during or before the fledgling 
period by selecting only birds ringed when fully grown to 
reduce age dependency of survival in the data. However, 
our estimated monthly survival probability of 0.89 
translates into an annual survival probability of 0.25. This 
estimate seems to be lower than or at the lower edge of 
annual survival estimates for the European robin given by 
other authors based on ringing and recovery data — Lack 
(1964): first year (including fledging period) 0.34, adults 
0.38, and migrants with unknown age 0.23; Fransson and 
Pettersson (2001): first year 0.21 and adults 0.35; and Siri- 
wardena et al. (1998): first year 0.40 and adults 0.42. 
Underestimation in survival probability can be caused by 



unaccounted heterogeneity in either survival or recovery 
probability (Carothers 1973; Gilbert 1973; Fletcher et al. 
2012). Such heterogeneity could be nonconstant survival 
probability over the year. For example, in the black- 
throated blue warbler, survival was reduced during the 
migration periods compared with the rest of the year 
(Sillett & Holmes 2002). 

In order to assess the effect of such nonconstant sur- 
vival in the data on our model estimates, we simulated 
data based on the assumption of lower survival probabil- 
ity during the migration period than during winter and 
summer. Fitting the model to these data showed that 
nonconstant survival primarily affects the estimates for 
recovery probability, but has no effect on bias or preci- 
sion of the distribution parameter estimates (Data S5). 
These simulations also show that the model in its present 
form is primarily designed for estimating bird distribu- 
tions rather than analyzing survival. 



726 



© 2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



F. Korner-Nievergelt ef al. 



Bird Distribution Model for Ring Recovery Data 



■ Central Europe □ Southern Europe 
□ Northern Africa 




a) 

"O 

S 0.6- 

o 
c 
o 

1 0.4- 

o 

Q. 



0.2- 



o.o 1 1 1 

Fennoscandia central Europe 

Figure 3. Mean of the estimated winter distribution (December- 
February) of birds ringed in Fennoscandia during all months with a 
sufficient high number of ringed birds (April, May, September, and 
October) and central Europe during the months May-August. The 
colors in the bars give the distribution of the birds among the three 
regions. The vertical lines give the 95% credible intervals. 

A further point that deserves discussion is that it is not 
possible to identify the sex from plumage in the European 
robin. As a consequence, we cannot include sex-specific 
migration in our model. However, it is well known that a 
higher proportion of the females are migratory compared 
with males. For example, at a Belgian study site, 100% of 
the females were migratory, whereas only 0-70% of the 
males were migratory (Adriaensen and Dhont 1990). 
Unaccounted heterogeneity in some parameters can 
sometimes cause bias in the parameter estimates in mark- 
recapture models (Fletcher et al. 2012). Therefore, we 
simulated data for males and females separately assuming 
different migration patterns for females and males. We 
then pooled the data (i.e., ignoring sex as in the real data) 
and analyzed the data with our model. The bias of the 
distribution parameters m,^ was, for the settings used in 
our simulation study, negligible (—0.003, if the mean of 
the sex-specific distribution parameters was assumed to 
be the true value), and the precision of the estimates was 
similar to data sets of equal sample sizes without distribu- 
tional differences between the sexes (Data S6). This means 
that the estimates corresponded, for our example, to the 
proportions of all birds (ignoring sex) in the different 
regions during the different seasons. It remains to be 
tested whether it is possible to account for such heteroge- 
neity by extending the model with a mixture model as 
proposed by Pledger (2000). 



While overall the model fitted the data well, the predic- 
tive model checking revealed some heterogeneities not 
accounted for in the model (Data S4). For example, pre- 
dicted numbers of reencounters in Fennoscandia during 
the summer were too high for birds ringed in Scandinavia 
in October. This may be because, during autumn, popula- 
tions from northeastern Fennoscandia are migrating 
through Europe and, thus, are ringed in Scandinavia and 
central Europe. These birds breed at places with much 
lower recovery probabilities than elsewhere in Fennoscan- 
dia. As a consequence, birds ringed during migration have 
a lower recovery probability during summer within 
Fennoscandia than birds ringed during the nonmigration 
period. As the model assumes that recovery probability 
does not depend on when the bird was ringed, we observe 
fewer recoveries than expected from those birds ringed 
during migration. That we do not see this discrepancy for 
Scandinavian birds ringed during spring migration may 
be because these northeastern breeding birds take more 
easterly routes in the spring. Such a loop migration for 
the northeastern breeding birds has been suggested by 
Pettersson et al. (1990) who analyzed morphological mea- 
surements of robins in southern Sweden and at different 
places in the wintering area. 

Literature that looks at the proportions of birds from 
different populations that stay in different regions over 
the year is sparse. Schifferli (1961) wrote that only 5-10% 
of the robins breeding in Switzerland stay there over win- 
ter. Among 46 individual offspring of a population from 
southern Germany, 10 (22%) did not show migratory 
restlessness (Biebach 1983). Our model estimated this 
proportion as 36% (12-63%) for those birds ringed dur- 
ing May and 19% (0.8-48%) for those birds ringed dur- 
ing June and July. Our model estimates are slightly higher 
than those in the literature. In Belgium, all females left 
the study plots (of 10-70 ha size) in winter, whereas of 
the males only 70% left the woodland plots and 0% the 
park and garden plots. It is not known what proportion 
of these individuals left central Europe. However, the 
study shows that, especially in urban areas, robins may 
stay within less than 5 km of the breeding place to spend 
the winter. Because we discarded recoveries within less 
than 5 km distance from the place of ringing, our esti- 
mates of sedentary birds may be too low. Despite this 
underestimation, we have a higher proportion of seden- 
tary birds than earlier studies. This difference may be due 
to a long-term change in migration behavior or simply 
due to uncertainty in all of these estimates as the number 
of ringed birds during the breeding season is low. 

The proportions of birds wintering in northern Africa 
are consistently higher for those birds ringed in Scandina- 
via compared with those ringed in central Europe 
(Fig. 3). This means that many Scandinavian individuals 
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overtake central European ones during autumn migration. 
This may indicate a leap-frog migration at least in parts 
of the population (Salomonsen 1955). 

The motivation for this study was to develop a method 
to correct the ring recovery maps given in bird migration 
atlases for spatially heterogeneous ring recovery probabili- 
ties. The model we present here is an extension of previ- 
ous bird distribution models used by Bauthian et al. 
2007, Thorup and Conn 2009, and Korner-Nievergelt 
et al. 2012. We extended these models by including a sea- 
sonal dependence of the bird distribution without includ- 
ing a year dependence of any model parameter. The 
independence from year makes our model suitable for 
analyzing data sets that have been collected over a long 
time period (>50 years) with few observations per year. 

The model enables quantification of the temporal and 
spatial distribution of bird populations as well as migra- 
tory connectivity. Based on the experience from this 
study, we recommend including as many individuals as 
possible in the data, including data from different species. 
As the assumption that recovery probabilities are equal 
between different species may not be realistic, species-spe- 
cific recovery probabilities may need to be estimated with 
partial pooling between the species. In addition, other co- 
variates for recovery probability, for example human den- 
sity (Korner-Nievergelt et al. 2010a), may be added to the 
model. Finally, the regions should be chosen small 
enough so that it can be assumed that ring recovery prob- 
ability is homogeneous within each region. 

The application of this model to combined data of a 
large number of species with a high number of regions 
requires that the computation is fast. The models pre- 
sented here needed around three minutes to be fitted to 
the present data set on an intel COREi5vPro notebook. 
This is fast enough to be confident that fitting a similar 
model to data of more regions and more sets of birds 
released will be feasible. 

Finally, we would like to note that the Bayesian frame- 
work allows the inclusion of information from other data 
sources, such as modern tracking methods, either as prior 
information or as an integrated model, similar to the 
integrated population models (Abadi et al. 2010). Data 
sets from different sources often differ in quality, for 
example ring recovery data contain limited information 
about many individuals, primarily after death, whereas 
tracking data contain more information of only a few 
individuals that survived until the data could be down- 
loaded. The combination of such different data sources 
has the potential to give a more comprehensive insight 
into the migration behavior of the species than two sepa- 
rate analyses. Our study contributes to the development 
of formal methods for the combination of different data 
sources that has been recognized as being important for 



F. Korner-Nievergelt et al. 

the study of migratory connectivity (Norris et al. 2006; 
Baillie et al. 2009; Fiedler 2009). 
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Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Data SI. R and JAGS code for fitting the distribution 
model to ring reencounter data. 

Data S2. The aim of this simulation study was to assess 
the sample size required to estimate the proportion of 
birds of the different sets of birds (if) being in the differ- 
ent regions k during each season q, mijkq. 
Data S3. This document contains two figures that show 
the 95% credible intervals of all proportions presented in 
the barplot of Fig. 2 of the main text. 
Data S4. Diagnostic figures from posterior predictive 
model checking for two models: (1) model assuming con- 
stant recovery probability over time (not presented in the 
main text), (2) model assuming season-dependent recov- 
ery probability (as presented in the main text). 
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Data S5. The aim of this simulation study was to assess 
the influence of nonconstant survival on the estimates of 
the model parameter, as our model assumes constant sur- 
vival probability. 

Data S6. The aim of this supplementary material is to 
assess the effect of sex-specific migration on the model 
estimates when sex is not identified. 

Fig. SI. Mean and 95% interval of the posterior predic- 
tive distributions of the number of recoveries [log 
(R new i j kq +l)] centered around the observed number of 
recoveries in each region for each season (bottom axis) 



for birds ringed in Fennoscandia during a specific month 
(left outer axis). 

Fig. S2. Legend as in Fig. SI but for birds ringed in cen- 
tral Europe. 

Fig. S3. (A) Estimated proportions of birds ringed in Fen- 
noscandia during different months (indicated on the left 
outer axis) being in the four different regions (indicated 
on the bottom) during 8 different seasons (x-axis). (B) 
Legend as in Fig. S3a but for birds released in central 
Europe. 
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